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This article considers a quaternary direct-sequence code-division multiple-access 
(DS-CDMA) communication system with asymmetric quadrature phase-shift- 
keying (AQPSK) modulation for unequal error protection (UEP) capability. Both 
time synchronous and asynchronous cases are investigated. An expression for the 
probability distribution of the multiple-access interference is derived. The exact 
bit-error performance and the approximate performance using a Gaussian approxi- 
mation and random signature sequences are evaluated by extending the techniques 
used for uniform quadrature phase-shift-keying (QPSK) and binary phase-shift- 
keying (BPSK) DS-CDMA systems. Finally, a general system model with unequal 
user power and the near-far problem is considered and analyzed. The results show 
that, for a system with UEP capability, the less protected data bits are more sen- 
sitive to the near-far effect that occurs in a multiple-access environment than are 
the more protected bits. 


I. Introduction 

In the design of a wireless communication system, feedback between the transmitter and receiver 
regarding the channel condition is useful for adapting the radio transmission rate to match the channel 
conditions [1,3-5,19]. When the channel condition is good, the data rate is increased, while when the 
channel condition is bad, the data rate is decreased. However, in some cases the transmitter does not 
know the condition of the channel and still desires to match the data rate to the channel. In this case, 
modulation and demodulation techniques are needed that allow more data to be transmitted when the 
channel is good and less when the channel is bad, without the transmitter knowing in advance the 
condition of the channel. 
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Consider for example the transmission of an image. Suppose that there are two modes of operation 
at the receiver with respect to high and low signal-to-noise ratios (SNRs). The two modes have different 
demodulation and decoding strategies according to two different rates and image qualities. In the high 
SNR mode, the receiver can demodulate and decode the data at high rate (or full rate) and recover the 
image with its high quality. In the low SNR mode, the protection available with coding and modulation 
is not adequate to protect all the data. However, it may be possible to decode only a subset of the bits 
that have higher error protection. In this case, the receiver demodulates and decodes the data at a lower 
rate, and recovers the image with lower quality as compared to the high-quality image. 

In a wireless network, the channel condition can vary for several reasons. One reason is just the 
change in the distance between the transmitter and the receiver. Another reason is that the multiple- 
access interference produces time- varying channel conditions. 

The key idea in designing such a system is to introduce modulation and coding schemes that pro- 
vide different error protection to different classes of data. The earlier work on multicasting [6,11-16] 
and unequal error protection (UEP) [8,18,20] examined such a system in the case of a mobile network 
downlink. This idea is essential when different portions of the source do not contribute evenly to the 
overall quality of the decoded information. The UEP technique is a simple and efficient method to satisfy 
such a requirement. The basic idea is to use a constellation with non-uniformly spaced signal points in 
the modulation scheme. The non-uniform nature of such a constellation results in different distances be- 
tween sets of signals and provides different levels of reliability against noise and interference and, hence, 
unequal error protection on different bits of a symbol. An asymmetric quadrature phase-shift-keying 
(AQPSK) constellation can be regarded as the simplest modulation scheme to provide the system with 
UEP capability. 

In [9], a quaternary direct-sequence code-division multiple-access (DS-CDMA) system is analyzed and 
an expression for the SNR is determined. However, the exact bit-error rate (BER) performance is not 
derived. In [7] and [17], the case of binary DS-CDMA with random signature sequences is investigated 
for binary phase-shift-keying (BPSK). Also, the Gaussian approximation to the interference is used to 
approximate the performance. In this article, we derive the exact BER for a quaternary DS-CDMA system 
and also derive the approximate BER using a Gaussian approximation to the interference for AQPSK. 
We consider a direct-sequence spread-spectrum modulation technique with asymmetric QPSK modulation 
that allows higher data rate transmission if the channel is good and a lower transmission rate when the 
channel condition is poor. We analyze the performance of a quaternary DS-CDMA communication using 
AQPSK modulation over an additive white Gaussian noise (AWGN) channel, with a correlation receiver 
that is coherent to the desired user. We look at both the cases of specific and random signature sequences 
being used in the system. 

This article is organized as follows. In Section II, the system model is introduced. In Section III, we 
derive the exact BER performance of the system. This also includes the derivation of the probability 
density function (pdf) of the multiple-access interference (MAI). A numerical example is given to illustrate 
the performance using a specific set of signature sequences. In Section IV, the random signature-sequence 
case is considered. The Gaussian approximation is used to model the MAI, and the approximate BER 
performance is obtained. In Section V, we generalize the signal model and examine the near-far effect 
on the system performance. 


II. System Model 

In this section, we describe the mathematical model of an asymmetric QPSK modulation system 
and characterize the receiver output. We consider an extension of the model described in [9] for asyn- 
chronous quaternary DS-CDMA. The model is shown in Fig. 1. The difference from [9] is that we consider 
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•J2P cos(/?)a, , (f) cos(2 nf c t + O t ) 



Fig. 1. A quaternary DS-CDMA communication system model. 


asymmetric QPSK so that the in-phase (I)-channel and quadrature-phase (Q)-channel bits have unequal 
energy. Suppose there are K users in the system. The quaternary signal of the fcth user is given by 

Sk{t) = si(t) + sf?(t) 


where 


s k (t) = V2P ■ cos/3 • afe(£)&fe(i) cos(27r/ c t + 0 k ) 

( t ) = V2P ■ sin /3 ■ a® (t)b® ( t ) sin(27r/ c t + 0 k ) 

In the above expressions, P is the transmitted power, /3 is the angle of the signal points in the asymmetrical 
constellation, aj.it) and are the spreading signals of the I and Q channels, bj,(t) and b®{t) are the 

user information being transmitted in the I and Q channels, and 0 k is the initial phase of the fcth user and 
is assumed to be uniformly distributed over the interval [0 , 27r] . The modulation constellation is shown 
in Fig. 2. In this scheme, we choose 0 < (3 < 7t/4. 

The information being transmitted by user k is represented by 


b k(t)= b k,j ' Pr(t - jT) 

j = ~ oo 
oo 

b k(t)= b k,3 -Prit-jT) 

j =- oo 
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Fig. 2. Asymmetric QPSK constellation. 


where b T k ,■ , b k .- £ {±1}, T is the symbol duration, and 


Pr(t) 


1, 0 < t < T 
0, otherwise 


The spreading signals are expressed as 


a k(t)= a lj ■ ~ J T c) 

j = ~ oo 


H <*kj ' W ~ jT C ) 

j— -oo 


where a k j £ {±1} are the signature sequences for the I and Q channels, T c is the chip duration such 
that T = NT C , and ip(t) is the chip waveform which is nonzero for 0 < t <T C . In general, we can choose 
any pulse shape as the chip waveform. However, to simplify the analysis, in the following we will assume 
that the rectangular pulse is used as the chip waveform, i.e. , = pr c {t)- 

The receiver is assumed to consist of a simple correlator matched to the desired signal. We examine 
both the time synchronous and asynchronous cases. Even though the asynchronous case is the more 
realistic case of the two, the synchronous case is more easily analyzed than the asynchronous case. When 
considering channel coding using linear block codes, it is very difficult to analyze the asynchronous case 
due to the dependency of bit errors within one block. 

A. Asynchronous System 

We first consider the asynchronous case. The received signal is given by 
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K 

r(t) = ^2 s k(t- Tk) + n(t) 
k = 1 


K 

= \/2P cos/? • a£(t - - 7*) cos (27 r/ c (£ - r fc ) + 0 fc ) 

fc=l 


K 

+ ^ \/2Psin/? • - r k )b^(t - r k ) sin (27r f c (t - r k ) + + n(t) 

k=l 


K 

= ^^ cos P • - r fc) cos(2t xf c t + </> fc ) 

fc=l 


K 

+ ^ \/2P sin f3 • a^(£ - r k )b^(t - r k ) sin(27 rf c t + </>&) + n(t) 

k=l 


where n(t) is an additive white Gaussian noise with zero mean and two-sided power spectral density 
N 0 /2. The time delay of the fctli signal is represented by Tk and (j>k = Ok — 27r f c Tk (mod 2tt). 

The analysis here basically follows the methods in [9] and [2]. Consider the output of the correlation 
receiver for the first user. The output of the I-channel correlator for the data bit b[ 0 can be decomposed 
into terms corresponding to the desired signal, the interference, and noise as follows: 

Z[ = f r(t)a[ (t) cos(2nf c t)dt 
Jo 


V2Pcos,f3 ■ a[(t)b[(t) cos(27t f c t)a[(t) cos(2irf c t)dt 


+ f V / 2Psin/3 • af (t)bf(t) sin(2n f c t)a[(t) cos(2-irf c t)dt 
Jo 



V / 2Pcos/3 ■ a{(t — Tk)bl(t — Tk) cos(27t f c t + 4>k)a[(t) cos(27t f c t)dt 


+ 


K 

E 

k=2 


r T 

/ V2Psm/3- 
Jo 


a2 {t - T k )b^(t - T k ) sin(27 rf c t + 4>k)a[{t ) cos(27 if c t)dt 


n(t)a[(t) cos(2irf c t)dt 


In the above expression, the component due to the desired signal is 
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A= / V / 2Pcos/3 • b[(t) (a{(t)) 2 cos 2 (2nf c t)dt 
Jo 

r T 1 

= V2P cos f3-b{ 0 - [l + cos(47t/ c £)] dt 

Jo 2 

= T \fpj2 cos P ■ b[ 0 


where the double frequency term is negligible since we assume f c (T c ) 1 . Because of the assumption 
of coherent reception, the component of the I-channel correlator output due to the Q-channel signal is 


B = 



V2Psin/3 ■ af(t)bf(t) sm(2nf c t)a{(t) cos(27r f c t)dt 


= V2P sin p ■ b® Q 



af(t)a[(t) sin(27T f c t) cos( 27 t f c t)dt 


= 0 


The interference component of the I-channel correlator output is given by 


C k 



VzPcos /3 ■ a k {t — T k )b k (t — Tfc) cos(27r/ c t + (f>k)a{(t) cos(2irf c t)dt 


f T 

= V2Pcos/3 / b k (t — Tk)al(t — r k )a{(t) cos(2tt f c t + 4>k) cos{2tt f c t)dt 

Jo 


r 1 r 

V2P cos Pj b I k (t-Tk)a I k (t-Tk)a[(t)- cos(</> fc ) + cos ( 4 tt / c £ + p k ) 


dt 


rT 

= \/ P/2 cos P ■ cos(4> k ) / b k (t - T k )a k (t - T k )a{(t)dt 
Jo 


= \/Pj2 cos P ■ cos {<j> k ) \bi + 6fc >0 -Rfcfi(Tfc) 


(1) 


where the time cross-correlations Rj/pr) and RjJpr) are defined as [9] 2 


2 Note that the “hat” notation on the cross-correlation functions is used to denote the correlation over the complementary 
(with respect to the symbol duration) portion of the integration interval. 
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Letting 


R k,i( T ) = [ a k(t ~ T ) a i (t)dt 
Jo 


Rk,i( T )= / a[{t - r)a\{t)dt 


"(6^, t ^)-T 1 b I k _ 1 R I k I 1 (Tk) + bl^Rl ^Tk) cos </> 


then Eq. (1) can be written as 


C fc = T^P/2 cos p • /ft (b{, T k , (f> k ) 


The total I-channel interference is then 


22 C k = T sfPj2 cos P 22 Ha , n , 4>k) 


where b 1 = (6| _i, b T 2 0 , ■ ■ • , b r K _j. b J K 0 ). Similarly, the component of the I-channel correlator output due 
to a Q-channel interferer is 

[ T 

D k = V2P sin/3 • a^(t - r k )b^(t - r fc ) sin(27 xf c t + (j> k )a[{t) cos(27 rf c t)dt 
Jo 

( T 

= V2P sin/3 / b®(t- T k )a®(t - T k )a[{t ) sin(27r/ c t + <j> k ) cos(27 xf c t)dt 

Jo 

fT 1 

= V2P sin p / b^{t-T k )a^{t-T k )a[{t)- sin(^ fc ) + sin(47r/ c t + <f> k ) dt 


f T 

= \JP / 2 sin/3 • sin(0 fe ) / b^ (t - T k )a^ (t - T k )a[(t)dt 
Jo 


2~P]2 sin/3 • sinOfc) b* + b^ 0 R^[(T k ) 


where, analogously to Eqs. (2) through (4), we define 


i? fcJ( T )= f a 2{t-T)a\{t)dt 

Jo 


R ki 0) = / a fe 0 - r K 7 
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and 


C ( 6 ?,r, 0 ) =r- 1 [6g. 1 i2gl(^)+6g 0 ^(7 1b ) 


sin </> 


Then Eq. (5) can be written as 


D k = T y/P/2 sin/3 • 1^ (ft®, T kl 0 fe ) 


and hence the total Q-channel interference is 


D fc = T\J~Pj2 sin /3 53 ($, Tfc, 


k—2 


k—2 


where 6 ^ 0 > ' " -i^k o) ■ The n °i se component of the I-channel correlator output is 

n[ = f n(t)a[ (t) cos(27r/ c t)dt 

Jo 

Note that n{ is Gaussian with zero mean and variance NqT / 4. In summary, we have 

l~p K l~p K kp 

z{ = ry — cos /3 • + 53 T Y Y cos P ' 7 M (^fc» "He, </>fc) + 53 T y ^ sin ^ ' J M (^fc ’ Tfe ’ ^ fc ) + n i 


rV? { fo lo cos / 3 + cos / ? 53 i M i^ki T k,<t>k) + sill /3 53 33 (&fe,'Tfe,^fc)| +«1 

l k—2 k—2 ) 


= V ^ [ COS ^( 6 1,0 + / /)] 


where 


K 

h = 53 T K i (^> T *’ + tan /3 • 33 (fefc , T-fc. 0fc) 


k—2 


Similarly, for the Q-channel correlator output we have 


z? = 


{ K K 

b i , o sin /3 + cos /3534a {^k^kAk) + sin/353 7 3? {^k^k^kj 

/c— 2 fc=2 


Ty^ [sin/3 (6^ 0 + Jq) 
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where 


K 

J Q = J2 T k,l (&k> r k,0k) +cot/3/0 ($,T k , </>*?) 
k= 2 

li Q i (bi,T,(t^ = T- 1 fefe-i-Rfc^Tfe) + 6fc, 0 ^(r fc ) sin(— </>) 

4 Q ? = T" 1 [&£_!<? (r fe ) + 6g 0 ^g?(TTfc)] cos</> 

and 

R l Q i( T )= [ a{{t-T)af(t)dt 

R l Q i( T )= J a k(.t- T )a?(t)dt 

R kf( T )= [ a k( t ~ T ) a ?( t ) dt 

R kf 0) = J a k (t - t K q 

Also, is Gaussian with zero mean and variance NqT/A. 

B. Synchronous System 

For the synchronous case wherein Tk = 0 for all k = 1, 2, • • • , K, the received signal is given by 

K 

r(t) = 5>(t) + n(f) 

fc=l 

With arguments similar to those in the previous section, we have the correlation receiver outputs 

K 

z{ = n[ + Ts/pj2 cos/3 6( 0 + 'Yh I l,\( h k,Qk) 

k = 2 
if 

= n? + T/P/2 sin/3 ^ 7^(b fc , 0 fc ) 

fc=2 

where ?r{ and ny 1 are Gaussian with zero mean and variance NqT/ 4, = (6f. 0 ,b ® 0 ), and 
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Jfc,i( b fcA) = illiKo,Ok)+tznp.i%[(b% 0 ,ej 
4 Q i( bfc, 0fc) = K 0. 0*) + cot 0 ■ I l? (6g 0 , 0 fc ) 


with 

= T_1 ' kfe.o ' ftk?i( 0 ) ' cos(Ok) 
iZ (bi 0 , Ok) = T- 1 ■ b% 0 • i#'( 0) • sin(0fc) 

4 Q i (bio, 0 k ) = T - 1 • ft' 0 • i2g(0) • sin(— 0^) 

J??(&? 0 A) = T_1 ' &?o ' <?(«) ' cos(0 fc ) 

Since there is no delay between users, we have 

rT N-l „T C 

#m(°)= / ai(t)a[(t)dt = a lj a lj ip 2 (t)dt 

Jo j=0 Jo 

p J\[ Ji 

^fej(0)= [ al (t)a[(t)dt = a k,j a l: [ tp 2 (t)dt 

Jo j=Q Jo 

p 2 . p< 

R k Q i( Q )= [ a k(t)a? (t)dt = a lj a ?,j f 

Jo j=0 Jo 

/■T JV -1 r T c 

R k ®( °) = / a fc (t)af(t)dt = ^2 a^jCiZ / ip 2 {t)dt 

Jo J=0 Jo 

Furthermore, since we use a rectangular chip waveform, that is, = Pr c (t ) , then / Q Tc ip 2 (t)dt = T c , 
and we can further simplify the above expressions as 

JV-l 

<i(0) = 2c ^ a fc,j a i,j 

i= o 

JV-l 

^j(O) = T c Y, <$Aj 

3=0 
N - 1 

<?(o) = 

i=o 

N-l 

<?(0) = ^ £ °fcj 
0—0 
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III. Exact Performance Analysis 


Our goal is to analyze the bit-error rate (BER) of such a system. In order to find the exact BER, 
we need to find the probability distribution of the interference. In this section, we derive the pdf of the 
interference for both synchronous and asynchronous cases. 

A. Average Probability of Error 

The average probability of bit error is given by 

Pe=\(P I e+P?) 

where P[ and P® are the average probabilities of bit error of the I and Q channels, respectively, and are 
evaluated as follows: 

Pi = i{Pr (Z-[ < 0|6' 0 = +1) +Pr (Z[ > 0|6{ >o = -l) } 

= i {Pr (r^Pj2 cos 0(1 + I,) + n[ < o) + Pr (t^/P/2 cos 0(- 1 + I/) + n[ > o) } 

Pr ( = 2P < -1 - 1 + Pr 

\T-\fPj2 cos/3 ) 

= ^|Pr (rif + 1 1 < -1) + Pr (m + // > 1) | 

= ^|l - Pr (-1 < rn + h < 1) | 

where nj = n[/(Ty/P/ 2 cos 0) is Gaussian with zero mean and variance (2 /iVg) -1 , E l = PT cos 2 (3 = 
E s cos 2 (3 , and E s = PT. Similarly, we have 

P? = \{l-Pv(-l<n Q + I Q <l)) 

where uq = nf /(T \J P/2 sin 0) is Gaussian with zero mean and variance (2E^ /No) -1 , and E® = 
PT sin 2 0 = E s sin 2 (3. 

In order to evaluate Pj. and P0 , we use the characteristic function method in [2] to compute these 
probabilities. In order to compute P ei we need to know the probability distribution of the sum of the 
noise and interference. We first obtain the characteristic functions of the random variables, and then 
derive P f r and P0 from the characteristic functions. 
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Let < f>„ J (u), <!>/, (i’) , and $/(u) be the characteristic functions of n/, 7/, and I = m + //. Note that 
they are even functions (d>(u) = <£(— u)), and <I>j(i>) = by the independence of n/ and 7/. 

The probability needed to compute P/ is obtained as 

Pr (—1 < m + 7/ < 1) = J fj(x)dx 

= 2 [ fi(x)dx 
Jo 

-'Ml: £•**-’-*)* 


1 / /-OO 


$/(u) cos (vx)dv ) dx 


io \J o 


/ d>/(u) (/ cos (vx)dx'j dv 


$/(u)u i sin(w)(iu 


The characteristic function of the interference in the I-channel, <h/(u), can be written as 

®l(v) = 4>„ 7 (u)4>j 7 (v) = $m{v) - $ ni {v) + $ ni (u)$/ 7 (w) = 4> nj (v) - $ nj (v) [l - <&/ 7 (u)] 


thus, 


2 7*00 2 /‘OO 

Pr (— 1 < m + 7/ < 1) = — / ^ > nj(u)u _1 sin(u)du / sin(u)$ nj (u) [l — $/ 7 (u)l du 

71 Jo n Jo 


where 




$ n/ (u) =exp ( -H^v 2 




The average probability of error of the I-channel is then given by 


p e = \ ~ ^ Pr (- 1 < «/ + 7/ < 1) 


1 1 

2 7 r 


1*00 POO 

/ i> _1 sin(v)& ni (v)dv H — / u _1 sin(u)$ nj (u) [l — $/ 7 (u)] du 

do 77 Jo 


= Q\\HJ L 


2El 1 


iVo / 7T do 


v 1 sin(u)$„ J (v) [l — $/ 7 (i>) dv 
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where 


Q{x) 



x > 0 


Similarly, let <& n<3 (^), and < f>Q(ti) be the characteristic functions of rig and Iq. Then we have 


2 r oo 2 poo 

Pr (-1 < riQ + Iq < 1) = - / v sin(w)$„ Q (v)dv / f _1 sin(w)$„ Q (v) [l - $/ Q (u)] dv 

71 Jo 77 Jo 


and 


1 1 f°° 1 

P e Q = / u _1 sin(u)$„ Q (v)dv + - / v _1 sin(w)$ nQ (v) [l - $/ Q (u)] dv 

z 77 Jo 77 Jo 


= Q 



v 1 sin(u)$ nQ (u) [1 - $ Iq (u)] dv 


Therefore, the average probability of error is given by 



+ <5 



1 f°° 

+ ^] Q u _1 sin(?;) {$„,(«) [l 


$h 0)] + ®n Q (v) [1 


®i Q («)] } dv 


Note that by representing the error probability in this way it is clear what the contribution to error 
probability is from noise and interference. When there is no MAI, i.e., the single user case, the MAI 
term in the above expression is zero, and the probability of error is the same as in the case of an AWGN 
channel. In general, the MAI term in the above expression does not have a closed- form solution and needs 
to be evaluated numerically. However, in order to evaluate it numerically, we need to find expressions for 
®h(v) and ®i Q {v). 

B. Asynchronous Case 

Here we begin to derive the characteristic function of the interference in the asynchronous case. The 
I-channel interference is given by 


It = 


K 

(^fe> T fc’ fa) + tan p ■ T k , <l> k ) 


k = 2 


where 
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( 6 ) 


( hl,T k ,<j>k ) 


c °feOfc) b j 


T 


r M ($,Tk,<t>k) = S111 y ^ [^,-i^fc,i(^fc) + &fc,o^M( T fc) 


Now consider IT C < r k < (Z + 1)T C . In this case, we have 

Rk.iiTk) = Ci\{l - iV)i?^(r fc - ZT C ) + Ci\{l + 1 - 7V)i?v,(r fc - ZT C ) 
^fi(r fc ) = Ci^R^n - ZT C ) + Cft(Z + I)!?.,/, (rfc - IT C ) 


( 7 ) 


where R^(t) and R^(t) are the autocorrelation functions of the chip waveform defined as 


Ri!,{t) = / - r)dt 


R^(t) = [ ip(t)ip(t + T c - r)dt 
Jo 


Similarly, 


R^in) = ($1(1 - N)%( Tk - IT C ) + <$[(l + 1 - N)R^r k - IT C ) 


RZ in) = cZmArk - ITc) + cZ(l + 1 )R^(T k - IT C ) 


In the above expressions, C k \(l) and (^{(Z) are the aperiodic cross-correlation functions defined as 



' ^~\N — 1 — l 

2^=o 

a k,j a i,j+l ’ 

0 < l < N - 1 


II 

>-c “ 

o' 

JV — 1+/ 
2^3=0 

- 

1~N <l<0 

(8) 


,o, 


\1\>N 



' N — 1 — l 

3=0 

I Q 

a k,j a i,j+n 

0<1<N-1 


a 

H 

^~\N — 1+/ 
3=0 

I Q 

°k., i a 7..r 

1~N <l<0 

(9) 


,o, 


\1\>N 



Here {a{. ,■} and {a® ■} are the spreading sequences of the I and Q channels of the fcth user. The charac- 
teristic function of // is given by 
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$h{v) = E{ exp(jvli)} 


= E < exp 


' / K 

jv f 7 M (&fc> T fe> + tan P ■ T ki , Tfe, </>fc 


\k=2 



K 

n ^ {exp jv (jl*! (b_i, T k , <i>k) + tan/3 • r fc , X)) | 


k—2 


K 

n 

k—2 


2tt T \ 4 / V 4 . 


XX jf jf exp [jv [l k l (bl, t, cj>) + tan /3 • I%[ (b® , r,<j>)^ j 


^ £ 


drdcj) 


K i /.(i+l)T c 

n^EE/ £ / exp 

fc =2 h Q i =0 ^ ZT <= 

l -fe -fc 


COS cj) 

jv ^ r - 


K ,- i { C " i ( 1 - N ) R ^( t - IT c ) 


+ 0^(1 + 1 - N)Rj/;(t - IT C )) 


+b {, o ( C&m^r - IT C ) + C&(1 + 1 )i^(r - ZT C )) 


+ jv tan/3 ^ [&£_i (cg(Z - iV)^(r - IT C ) + C^(l + 1 - JV)ify(r - ZT C )) 


+6go (Cm( 0X(^ - /r c ) + Cg[(* + l)i^(r - ZT C )) 


drdcj) 


(10) 


With further simplification (see Appendix A), we obtain 


K ( N-l / 8 

$// 0 ) = II { X [' 52 f ( v ' i l ’ 9 i ( l ), h i { l),t 

k—2 l 1=0 \ i = 1 


where 


/*2"7t i* T c p ^ 

f(v;l,g(l),h(l),a)=—— / cos {- (cos</> • p(Z) + asin</> • h(l))R^(r) 

c J 0 « 0 


+ (cos</> • g(l + 1) + a sin <j> ■ h(l + 1 ))i?. i /,(r) | drdcj) 


and 
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to 

II 
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MO 

II 

<N 
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c. 
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MO 
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MO 

56(0 = 

C. 

MO 

to 

-A 

II 
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MO 

tO 
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II 
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If we consider a rectangular chip waveform, we 
(see Appendix A) 


0 ®[, a\ = tan fj 
Ci> «2 = -tan /3 
Ci> 03 = tan /3 
Q®\, 04 = — tan /3 

C, 05 = -tan /3 
Ci> 06 = tan /3 
Ci> 07 = — tan /3 
Og = tan /3 

l further simplify the function /(v; l,g(l), h(l), a) as 


f(v;l,g(l),h(l),a ) = J sinc {^^( cos 0 (g(Z + 1) - ff(0) + asin<£(/i(Z + 1) - *(0))} 

• cos{^(co8 0( 5 (i + 1) + </(0) + osin</)(/i(/ + 1) + /i(0)) | d(f> 


This expression is simple to evaluate numerically, which allows us to compute the characteristic function 
and the average error probability. From the above expression, we see that the characteristic function of 
the interference does not depend on the signal energy or SNR. The advantage is that we need to compute 
the characteristic function of the interference only once, and it can be applied to different SNR values to 
compute the probability of error. 

C. Synchronous Case 

For the synchronous case, the derivation is similar to the asynchronous case. The expressions for 
the bit-error probability for the I and Q channels are the same as for the asynchronous case. The only 
difference is in the expressions for the characteristic functions of the interference. These are given by 

$/M) = n \ T~ / cos C [ cos 4* ' R k*i(0) + tan/3 • sin cf> • -R/?{(°)] ) 

fc = 2 77 

+ cos ^ [cose/) • — tan/3 • sin^> • R® {(0)]^ d</)j 

K y" -j 27T 

$ /g(»)= n r cos (t [cos </>• R^(0)+ cot/3- sin </> • i?'p(0)]) 

T cos [ cos 0 * Rfej (0) - cot/3 • sin <j> • CC 0 )]) d(/>j 
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D. Numerical Examples 


Here we present a numerical example for the asynchronous case. In [2], the average error probability 
for a direct-sequence spread-spectrum multiple-access (DS-SSMA) system with symmetric QPSK modu- 
lation is investigated. The performance is evaluated using auto-optimal, least side-lobe energy (AO/LSE) 
sequences [10] as the spreading codes for the users in the system. For the quaternary system, the spread- 
ing factor is chosen to be N = 127, and there are 9 pairs of codes listed. In each pair of codes, the I- and 
Q-channel sequences are the reverse of each other. The AO/LSE codes for N = 127 are listed in Table 1. 

Each row represents a pair of codes. The generator polynomial coefficients are denoted by H and H 
in octal. The initial values in the shift registers are denoted by a 0 and 1 . The in-phase interference 
characteristic function from the second user to the first user using the above spreading codes with (3 = 7t/4 
is shown in Fig. 3. Since in the symmetric constellation the I- and Q-channel signals have the same power, 
the resulting characteristic functions of the I- and Q-channel interference are the same. Therefore, we 
show only the characteristic function of the I channel. For (3 = 7t/8, even though we use mutually reversed 
spreading codes for in-phase and quadrature-phase components, the characteristic functions are different. 
This is due to the unequal power of the I- and Q-channel signals in the asymmetric constellation and the 
cross-correlation nature of the spreading codes. The characteristic functions of the interference from the 
second user to the first user when (3 = tt/8 are shown in Fig. 4. The average probability of error when 
the number of users varies from 1 to 9 is shown in Fig. 5. The performance is worse than the symmetric 
case as shown in [2] . This is because the performance is dominated by the Q-channel performance, which 
is bad due to the low transmitted power. 


IV. Approximate Performance Analysis 

As seen in the previous section, the expressions for the interference are very complicated, and the 
evaluation for the exact performance is computationally tedious. Also, as in the numerical example, the 
results are for a specific set of signature sequences. One way to solve this problem is to use a Gaussian 
approximation to model the interference and to use random signature sequences in the analysis. Then a 
simple approximate expression for the BER can be obtained involving only the signal-to-interference-plus- 
noise ratio (SINR) and the Q function. In this section, we approximate the interference as a Gaussian 
random variable and assume random signature sequences. We find the variance of the interference and 
examine the approximate system performance. 


Table 1. AO/LSE codes (N= 127). 


H 

“o 

H~ 1 

V 

M 

L 

S 

211 

0010000 

221 

1001101 

17 

6 

2183 

217 

0000101 

361 

1111111 

15 

12 

2015 

235 

0001100 

271 

1000101 

17 

10 

2283 

247 

0010111 

345 

0110001 

17 

8 

2255 

277 

1110001 

375 

0101010 

19 

4 

2295 

357 

1110010 

367 

0110101 

17 

4 

2563 

323 

1110111 

313 

1000111 

17 

4 

2203 

203 

1101101 

301 

0010010 

17 

4 

2087 

325 

0000101 

253 

1101100 

19 

6 

2483 
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Fig. 3. In-phase interference characteristic function (A/= 127, (3 = k!A). 




Fig. 4. Interference characteristic functions (N = 1 27, (3 = 71/8): (a) in-phase and 

(b) quadrature. 
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Fig. 5. Probability of error for asymmetric QPSK DS-SSMA ( N = 127, p = jt/8). 


A. Asynchronous Case 

In order to find the approximate BER performance, we approximate the interference as a Gaussian 
random variable and find its variance. We first find the conditional variance of the interference, and then 
average over the random variables to find the variance. Therefore, we can obtain an expression for the 
SINR, and thus the approximate BER. 

The decision statistics at the output of the correlation receiver for user 1 are 

K 

Z[ = n[ + Ty/P / 2 • cos/3 • b[ 0 + y / P/2y^cos/3 • IT l 1 ■ cos(4> k ) + sin/3 ■ W® 1 ■ sin (<j> k ) 

fc = 2 
K 

= n i + T \J P/2 ■ cos /3 ■ bf 0 + y/P / 2 ^2 sin P • ■ cos (4>k) - cos /3 ■ Wl Q ■ sm((j> k ) 

fc = 2 

where 

Wl 1 = bl_ x ■ + bl i0 ■ 

W? 1 = b %_ r • R%\(r k ) + b % o • RZ(r k ) 

W? Q = b%_ , • <?(r fc ) + b% ■ R^(r k ) 

Wl Q = bl_, ■ Ri% k ) + b{ fi ■ Ri% k ) 
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To find the variance of the multiple access interference (MAI) of Z[ and Zf, we start by writing Z( in 
the form 


Z{ = n[ + T^P/2 ■ b[ 0 • cos/3 + IT 


where 


K 


K 


W = \/P/2 ■ cos / 3 ^ Wl 1 ■ cos(4>k ) + \JP /2 • sin /3 W® 1 ■ sin (fa) = W 1 + W® 


k—2 


fc = 2 


with 


K 

W 1 = \Zp/ 2 • cos p ^2 w k • cos(^fc) 

k—2 


K 

W Q = a/ P/ 2 • sin/3^ • sin(0/ c ) 

k—2 


The variances of W 1 and W ® are given by (see Appendix B) 

y„[w-] = 

Hence, the variance of the MAI in Z{ is given by 

Var[W] = Var[IT 7 ] + Var[IT Q ] = ~ l )^ PT c 

The SINR of Z[ is then 


SINK/ 


T 2 P/ 2 • cos 2 /3 
N 0 T | (AT-l)fVPT 2 

"T~ + 6 


6i3 s • cos 2 (3 


3N 0 + 2E S 


(K - 1) 
N 


12 E b ■ cos 2 /3 


3 -/V 0 + AE b 


{K - 1 ) 
N 
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where E s = NPT C is the symbol energy, and E b = (1/2 )E S is the average bit energy. Similarly, for the 
Q channel, the SINR of Z ® is given by 


SINRq 


6E a ■ sin 2 (3 


3Nq + 2 E s 


(K - 1 ) 
N 


12 E b • sin 2 (3 


3 N 0 + 4 E b 


{K - 1 ) 

N 


Then, the approximate BER can be expressed as 


pi 

r e,GA 



\ 

12 E b • cos 2 (3 

\\ 

(K - 1) 

3N Q + m K N J 


E 


Q 


e,GA 



\ 

12 E b ■ sin 2 / 3 

A 

3N 0 + AE b /m^ ^ 


B. Synchronous Case 

The analysis for the synchronous case is similar to that for the asynchronous case presented in the 
previous subsection. We can rewrite the decision statistic as 

Z[ = n[ + T \f~Pj2 - cos (3- b[ 0 + W 


where 


K 

w = Y^ w l 

k—2 

and 


= s/m (cos/3 • b[ fi ■ Rk,m ■ cos(6»fc) + sin/3 • b% 0 ■ R%[( 0) • sin(6> fc )) 

We want to find the variance of the MAI W l . Note that and R®[( 0) are both functions of {a{ •}. 

Thus, the variance of Wj. conditioned on {ufj} and 6k is 
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Var [W k | {a{j}, 9 k ] = Var sJP/2 ■ cos f3 ■ b{ 0 • #^(0) • cos 6 k \ {a{j}, 9 k 


+ Var 


\J~Pj2 • sin/3 • b® Q • i?^(0) -sin6» fc | {a[j},0 k 


P 

2 " 


(cos 2 (3 ■ cos 2 9 k ■ Var b k 0 ■ R 1 ^ i(0) | {a{j} 


+ sin 2 (3 • sin 2 9 k ■ Var b^ 0 ■ R k ^( 0) | {a{j} ) 


Because we assume random signature sequences, given {a[ •}, /^(O) and R k [{ 0) are independent 
identically distributed (i.i.d.) with pdf 


N 


p R (rT c ) = r + jsf 2 


— N 


for r = —N, —N + 2, • • • , iV — 2, A/”. Since b k 0 and Rj/i have zero mean and are independent, we have 


Var [b{ o • Rl 1 ^ 0) | {<,•}] = E (6{ %0 • i?"(0)) 2 I {<,-} 


= £ 


(&io) 2 E (i*"( 0)) 2 | {a{j} 


= 1 • ATT 2 


= iVT 2 


Note that even though Rj, 1 ^) depends on {a{ • }, the mean and variance do not depend on the particular 
realization of {a [ }. This is different from the asynchronous case. However, this property helps reduce 
the complexity of the analysis. Similarly, we have Var[&j? 0 • R® i(0)|{a( J }] = NT%. Therefore, the 
conditional variance of is 


Var \W k | 9 k ,{a[ tj }\ 


NPT 'l 
2 


(cos 2 (3 • cos 2 9 k + sin 2 (3 • sin 2 9 k ) 


Let 0 = (0i, • • • , 9k)- The conditional variance of W is then given by 


Var[fV | ©,{<,}] 


K 


E Var [Wl | 9 k ,{aij}] 

k = 2 


= (K- 1) 


NPT 'l 
2 


(cos 2 (3 • cos 2 9 k + sin 2 (3 ■ sin 2 9 k ) 
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Note that the above expression now depends only on 9 k - By averaging over 9k, the variance of W is 


I< 


Var [W] = £> fc [Var [Wl\6 k ]} 


k—2 


= (K~ 1 )Eg k 


NPT? 


(cos 2 /3 ■ cos 2 9k + sin 2 /3 ■ sin 2 9k) 


(I< - 1)NPT 2 


(cos 2 (3 • Eg k [cos 2 9 k ] + sin 2 P ■ Eg k [sin 2 9 k }) 


K - 1 )NPT 2 ( 2 a 1 .2 a 1 

v 2 c ' cos 2 /3 h sill 2 p ■ - 

2 2 


(. K - 1 )NPT 2 


Therefore, the SINR is 


SINR/ 


T 2 P/ 2 • cos 2 P 


N 0 T / 4 + (A' - l)7VPT 2 /4 


2NT C E S cos 2 /3 
N 0 NT c + {K- l)T c E s 


2E S cos 2 P 
iV 0 + ^ 


4Af, cos 2 P 


By approximating the MAI as Gaussian with variance (AT — l)AP/4, the approximate I-channel average 
probability of bit error is 


PLga = Q (\/s!nr7) = Q 


( 


4 Eg cos 2 P 

2(K-1) 
Ao + ^ 


Similarly, it can be shown that, for the Q channel, the approximate average probability of bit error is 
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Pe,GA ~ Q 


4 E b sin 2 (3 
, 2(X-1) 


V. A Generalized Model and the Near-Far Problem 

In this section, we consider a general model for the AQPSK DS-CDMA system. The main difference 
from the model in the previous sections is that the users can have different transmission power. This 
causes what is referred to as “the near-far problem.” We are interested in the near-far effect on system 
performance. 

A. Analysis 

In the general model, the in-phase and quadrature components are given by 

s{(t) = A k cos/3-a{(t)bl(t)cos(2irf c t + d k ) 


s k (t) = A k sin f3 ■ a® (t)b% (t) sm(2n f c t + 9 k ) 


where Ai,A 2 ,--- ,Ak can be different. Without loss of generality, let user 1 be the desired user. The 
correlation receiver output of the in-phase and quadrature-phase channels are 


z{ = \a x T ib[ t0 cos P + cos p^2^ll A (bi,T k ,(l) k ) + ,T k ,<f> k ) \ + n[ 


z i = \ A i T l b i,o sin/? + cos/3^^/^(&£,T fc ,0k) +sin/?^^/^ 1 (6^,r fc ,^fc)| + : 
l k=2 k=2 ' ) 


The average probability of error is given by 


P e = 2 (Pi + P?) 


p ‘=W- p {- i <A +, p i 


P ' Q “ A 1 ^ p I -1 < a §^ +1 ' £ 1 


where 
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9k,i( l ) = Ci> h k,i( l ) = Ok,i> a i = tan 0 
9 k , 2 (0 = d k,i , h[,2 (0=^i. «2 = - tan 61 

9 I k,^) = d k I n h{ t 3 (l) =e^ t [, a| = tan 0 

9k,S) = e Z l. a| = -tan6» 

9k, 5 (0 = 9k,: l . ftfc.s (0 = - a l = - tan 0 

5fc,e(0 = <?m> ft fc,6(0 = a l = tan (9 

9k,7( l ) = 9k, i, K,r{ 1 ) =9k,v a 7 = - tan 0 

Sk,8 (0 = tflfi, /ifc,s(0 = C’ a 8 = tan 61 
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Similarly, we have 


where 


and 


K 

$/«(«) = n <bjQ 0) 

1 k , 1 

fc= 2 


(v) = 

1 

AT — 1 8 , 

^ EE /(^ 

i9k,i( l ), h k, 


I «- 

^ 1^ 


1=0 

i = 1 x 




J- / 

9k, i( 0 

_ nQQ 
— u k, 1 ’ 

*gi(D = 

e IQ 

Q 

a l 

= — cot 

9 

S ? 2 (0 

nQQ 

~ u k, 1 ’ 

< 2 (0 = 

e IQ 

°k, 1’ 

Q 

= cot 9 


9k, 2,(1) 

nQQ 

~ u k,l ’ 

h Q k ,S) = 

e IQ 

u k,l’ 

Q 

a 3 

= — cot 

9 

9k, 4 (0 

_ nQQ 
~ u k, 1 ’ 

hlS) = 

e IQ 

u k,l> 

Q 

o^t = 

cot 9 


9k, S) 

_ nQQ 
~ u k,l ’ 

< 5 (0 = 

0 IQ 

u k,l> 

Q 

= cot 9 


9k, S) 

_ nQQ 
~ u k, 1 ’ 

h Q k ,S) = 

e IQ 

u k,l> 

Q 

a 6 

= — cot 

9 

9k,r( l ) 

_ nQQ 
~ u k, 1 ’ 

h Q k , 7 (l) = 

e IQ 

u k,l' 

Q 

aj 

= cot 9 


00 

_ nQQ 
~ u k, 1 ’ 
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B. Numerical Examples 

Here we show some numerical examples. To see the near-far effect on the error probability, we consider 
the cases where there are five users in the system and the desired user’s power is four times the power of 
the interferers, while the interferers have the same power — that is, P\ = 4 P 2 = 4 P 3 = 4Pt = 4 P 5 . Here 
the total interference power is the same as Pi. We compare it with the case when there are two users 
having the same power, i.e., Pi = P 2 . In this case also, the total interference power is Pi. 

Figure 6 shows the average probability of error for both the I and Q channels with (3 = 7t/8 and 
N = 127. Due to the unequal error protection for the I and Q channels by the modulation scheme, we 
can see that the I channel has much lower error probability than that of the Q channel. 

Figures 7 and 8 show the average error probability for the I and Q channels for the two cases when 
9 = 7t/8 and N = 127. As can be seen, in the case with the near-far effect, the performance is better as 
SNR increases. This is because even though the total interference power is the same, the effect of each 
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K= 9 
K= 1 


K= 9 
K = 8 
K = 7 
K = 6 
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K = 3 
K= 2 


K = 1 
AWGN 


Fig. 6. Probability of error for asymmetric QPSK DS-CDMA (A/ = 127, (1 = n/8). 



£ b //V 0 , dB 

Fig. 7. Probability of error for asymmetric QPSK DS-CDMA l-channel 

(A/ = 127, p = ji/8). 
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E b IN 0 , dB 

Fig. 8. Probability of error for asymmetric QPSK DS-CDMA Q-channel 

(N= 127, (5 = ti/8). 


interferer on the desired user is not the same due to the different correlation relations of the spreading 
codes. In this case, the interference effect is not four times that of any one interferer since it is unlikely 
that all interferers’ spreading codes have simultaneously large correlation with the desired user. 


VI. Conclusions 

In this article, the exact and an approximate BER performance were derived for a quaternary asym- 
metric QPSK DS-CDMA system. The variance and pdf of the MAI were analyzed. The results showed 
that the AQPSK scheme can provide a significant difference in the amount of error protection for differ- 
ent bits of a symbol. Therefore, it is advantageous to use AQPSK when designing a UEP system for its 
simplicity and efficiency. We also examined the near-far problem by generalizing the system model to 
the case where users have different transmit power. The results showed that the Q-channel (less power) 
is more sensitive to the near-far effect than the I-channel is in a multiple-access environment. 
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Appendix A 

Characteristic Function of I T 

To further simplify Eq. (10), we have 

K I 1 JV_1 C 2 ™ f T Till 

= EE ex p 3 v t (A) dTd(t> \ 

fc = 2 ( *=° ° ° k Q 1 J J 

where 

A = [cos -N) + b[ 0 Ci\{l)) 

+ tan (3 • sin - N) + &g o Cg[(0)] ■ R^{r) 

+ cos<j!)(6^ i _ 1 C'^i(/ + 1 — N) + blgCl^l + 1)) 

+ tan (3 • sin ^(6g_ x Cgf(2 + 1 - N) + &g 0 C$(i + 1))] • R 4 ,{t) 

To evaluate the summations over bj. and b k , we note that there are 16 cases for ( b{,b = (b ! k _ 1 ,bj, 0 , b 
b k0 ) as in Table A-l. 

Introducing the periodic cross-correlation functions 

Ci(0 = O) + O - N ) 

W,s) = O) - O - n) 
o%!s) = tf!!s) + tf!!s- N ) 

we evaluate A for the above 16 cases with the following results. 
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j ^<0 



Table A-1. Sixteen cases for (b*, b ®). 


Case 


b l, 0 

b Q 

°k,-l 

b Q 

°k, 0 

1 

1 

1 

1 

1 

2 

1 

1 

1 

-1 

3 

1 

1 

-1 

1 

4 

1 

1 

-1 

-1 

5 

1 

-1 

1 

1 

6 

1 

-1 

1 

-1 

7 

1 

-1 

-1 

1 

8 

1 

-1 

-1 

-1 

9 

-1 

1 

1 

1 

10 

-1 

1 

1 

-1 

11 

-1 

1 

-1 

1 

12 

-1 

1 

-1 

-1 

13 

-1 

-1 

1 

1 

14 

-1 

-1 

1 

-1 

15 

-1 

-1 

-1 

1 

16 

-1 

-1 

-1 

-1 


Case 1: 

A x = cos<j> ■ 0^(1) + tan/? • sin</> • R^(t) 

+ cos(j) ■ 9 I k I 1 {l + 1) + tan/3 • sin<p • 9 k I 1 (l + 1) R^(t) 

Case 2: 

A 2 = cos <j) ■ Oj. 1 ^ ( l ) — tan f3 • sin cj) ■ 9^[ ( l ) (r) 

+ cos</> • 9 I k I 1 {l + 1) — tan/3 • sin^ • 9^\{l + 1) R^(t) 

Case 3: 

A 3 = cos (f> ■ 9 ( l ) + tan /3 • sin </> ■ 9^[ ( l ) R^ (r) 

+ cos</> • 9 I kl {l + 1) + tan/3 • sinc/i • 9 k [(l + 1) R^{t) 
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Case 4: 


A 4 = cos (j) ■ 0 kl (1) — tan (3 ■ sin (j) ■ 0 kl ( l ) (r) 


Case 5: 


Case 6: 


Case 7: 


Case 8: 


cos</> • + 1) — tan/3 • sin <j> ■ 0 k [(l + 1) R^(t) 


^ cos <t> • 61" (0 - tan/3 • sin <f> ■ 0 k \(l) R^(t) 


+ cos (j) ■ 01/(1 + 1) — tan f3 ■ sin <j> ■ 0 kl (l + 1) R^(t 


^ cos</> • Ol/(l) + tan/3 • sin</> • 0^[(l) R^(t) 


+ cos (j) ■ 0j,/(l + 1) + tan f3 ■ sin <j> ■ 0 kl (l + 1) R^(t 


^ cos</> • Ol/(l) — tan/3 • sin</> • 0^/(1) R^(t) 


+ cos (j) ■ 0 kl (l + 1) — tan f3 ■ sin <j> ■ 0 kl (l + 1) R^(t 


^ cos <j> ■ (l) + tan /3 • sin <f> ■ 0 k [ (l) R^(t) 


+ cos (j) ■ 0 k l (l + 1) + tan f3 ■ sin <j> ■ 0 kl (l + 1) R^(t 


The latter 8 cases are simply the negative of the first 8 and thus we have 8 pairs of cases. For case 1 
and case 16, we have 

ex P + exp = e JV< ' Al/T) + e ~ 3V ^ 1,T) = 2 cos (^ 1 ) 

with similar results for the other pairs of cases. Therefore, we have 


k ( 1 JV-1 / 8 


fe = 2 t i=0 \i=l 
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If we further define 


f{v;l,g(l),h(l),a) = 




cos0 • g(l) + 


a sin 


KI))%{t) 


+ ( cos (j) ■ g(l + 1) + a sin (j> ■ h(l + 1 ))R^(t) 


drdcj) 


then the characteristic function can be written as 


K 


N-l 


$ 


i A 


( v ) = II 1 87v ( '52f( v ’ l ’9iQ),hi(l),a i 


k = 2 


1 = 0 


where 


0i(O = 0"iMl) = C’ Q i = tan/3 

52(0 = 0 I f / 1 ,h 2 (l) = ^,a 2 = - tan (3 
93(1) = ^fi,/i3(0 = ^,03 = tan/3 
54(0 = = ^.«4 = - tan/3 

55(0 = 0{\,h 5 (l) = 0^[,a 5 = - tan/3 
06(0 = 0fcfi,MO = ^J,ae = tan/3 
97(1) = ei^hril) = d^[,a 7 = - tan/3 
0s (0 = ^fcfi) ft 8(0 = = tan/3 


As can be seen, in order to evaluate the characteristic function, we need to evaluate f(v,l,g(l),h(l),a), 
which involves the computation of double integrals that can be complicated. We can further simplify this 
by integrating over r when considering the chip waveform ^(f) to be the rectangular pulse. In this case, 
R^(t) = T c — t and R^(t) = r. Now the integrand can be written as 
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E= cosj^ (cos (j) ■ g(l) + asin</> • /i(Z))i?^(r) + (cos <j> ■ g(l + 1) + asin</> • h(l + l))i?^,(r) j 


cos |^ ( cos ^ • g(}) + asin<^ • h(l))(T c — r) + ( cosc/> • g{l + 1) + asin^> • h(l + l))r j 


cos cos (j>((g(l + 1) - g{l))r + g{l)T^j + asin cj)^(h(l + 1) - h(l))r + h(l)T^j j 


cos (^cos 4>(g(l + 1) — g(l]) + asin cf)(h(l + 1) — /i(/))^r + ^cos 4>g(l) + asin 4>h(l)^T c | 


= cos ( Ft + G ) 


where 


F = - (cos cj>({g(l + 1) - g{l)) + as\\uj)((h(l + 1) - h(l))) 


and 


G= — ( cos 4>g(l) + asin cj)h{l))T c 


Next, we integrate E over r to obtain 


yJ cos (Ft + G ) dr = sin (Ft + G) 


1 

'FT r 


| sin ( FT C + G) — sin ( G ) j 


i^{ 2 sm (irrd cos Gft + g 


= sine ( Y FT ^) C0s i\ FTc + G 


where 


sine (a:) = 


sin(7ra;) 


TTX 


1 1 V 

— FT C = 

2tt 2tt T 


^cos (j)(g{l + 1) — g(l )) + asin <j)(h(l + 1) — /i(Z))^T c 


2t tN 


^cos 4>(g(l + 1) — g(l)) + asin cj>(h(l + 1) 
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and 


\ FT c + G= ^(cos <j>(g(l + 1) - g{l)) + a sin + 1) 


— (cos(j)hi(l) + a sin <^/i 2 (0)^c 


^(cos (j>(g{l + 1) +s(0) + a sin <j>(h(l + 1) 


Therefore, we have 


f(v;l,g(l),h(l),a ) = ^ sinc {^^ (cos<l)(g(l + 1) - g(l)) + asm(j)(h(l + 1) 


• cos ( cos4>(g(l + 1) + g(l)) + asin <f>(h(l + 1) + /i(/)))| d<f> 


Appendix B 

Variance of the Interference 

The derivation here follows the work in [7]. According to Eqs. (6) and (7), we can expand Eq. (11) as 


W? = [&£,_! • C{\{ lk — n) + bl 0 • c &( 7fc )] Ms k ) 

+ [bl - 1 • Cfcfi ('Yk + 1 - N) + bl o • c"( 7fc + 1)] i^(Sfc) 


(B-l) 


where S k = T k — ^ k T c and j k = [Tfe/T c J. If we use Eq. (8) to expand Eq. (B-l), we obtain 


w[ r = 


1k~ 1 

^ ' bk,-l a k,j-'y k +N a 
3=0 


N-l 




O a k,j-j k a l,j 


R 4 ,(S k ) 


+ 


Ik 


N-l 


b k ,-ia k ,i-'y k -i+N a i,i + b ko a k j_ lk _ 1 a 1 j 


3=0 


i=7fe+i 


i^(-sfc) 


(B-2) 


which can be further expanded to obtain 
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w. 


II 


7k — 1 N—2 

^ — l^k,j— 7^+JV^l,j "b y bk,0®k,j— 7fc ®l,j 

3=0 j=7k 


bl,O a k,N-* k -l a l,N-l 


^(Sfc) 


7fc — 1 


^fc, — l^fc,W— 7 ^ — 1 ^ 1,0 "P ^ , ^k y — l®k,j— 7 fc+AT^l ) j-(-l T 


3=0 


N—2 

E 

3=lk 


K,o 


fe>3-7fe°l,3 + l 


i^(Sfc) (B-3) 


Finally, the terms in Eq. (B-3) can be rearranged to obtain 

7k — 1 

W 7 *, = b k t —i y ' dkj-jk+N (^ a i,j-^4’(^k) + oi j j-|_ii?^,(S'fc)^ 

3=0 

N—2 

+ 4o E a fe,3-7k (< jR'lp^Sk) + + (*5*^)^ 

j=lk 


+ 4o a fc,Ar-7k-i a i,A r -i-^V’(‘^fc) + bj, _ 1 a I k N _^ k _ 1 a[ 0 R^,(Sk) 


(B-4) 


In order to reduce the complexity of evaluating Eq. (B-4), we consider it conditioned on the signature 
sequence of the first user {a{ ■} and the random variable 7*,, which is uniformly distributed on the set 
{0, • • • , IV — 1}. We condition on 7*, = % and {a{j} = {af^-}, and define a set of N + 1 random 
variables Slj, 0 < j < N, by 


Q, 


,-l a fc,3-7k+iV a l,3> 

K, o a hj-% & ij> 

< 

b k ,O a k,N-*/ k -l b l,N-li 
, b k -l a k,N-' y fc -l®l,0> 


j = 0, • • • , 7fc - 1 
j = Ik,-", N - 2 

j = N-l 
j = N 


Then Eq. (B-4) can be simplified to 


W, 77 = 


N—2 

E 

3=0 


^3 


Rip(Sk) + b{jb{,j+iRi[’(S k ) + flAr-ii? i /,(S'fe) + Hat I?,/, (-Sfc) 


(B-5) 


where the random variables Qj , 0 < j < TV, are mutually independent and satisfy Pr(f2j = +1) = 
Pr(flj = —1) = 1/2. If we further define /(s) = R^(s) + R^is), g(s) = R^(s) — Ri/j(s), the set Ti to be 
the set of all nonnegative integers i less than N — 1 such that a 7 jfi 7 i+1 = 1 and the set T2 to be the set 
of all nonnegative integers i less than TV — 1 such that a[ ^a 7 i+1 = —1, then Eq. (B-5) can be written as 


Wl 1 = Y, to,f{Sk) + E n i9(Sk) + n N -i R^(Sk) + n N R^(s k ) 

36 r.i 3'er 2 


If we let X]J = Y 


jer 


% Yl 1 


= Yje r 2 ^3 > ^fc 7 = fljv-i) and A jj = fljv, then we have 
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Wi 1 = n I'R^Sk) + A I'R^Sk) + X 77 /(S fc ) + Y^g(S k ) 

Similarly, W k X can be written as 

wf 7 = n £%(S fc ) + A FMSh) + X^f(S k ) + Y k QI g(S k ) 

with n^ 7 , A^ 7 , AT]® 7 , and Y^ r defined in a similar way. 

At this point, in order to simplify the notation, we ignore the superscript of W k , II*,, A k , X k , and Y k . 
The random variables 11^ and A k are uniform on {—1,1}, and X k and Y k have pdfs 

Px k (i) = + L , i G {— L, — L + 2, • • • , L — 2, L} 

( M \ 

PY k (i)= \i + M 2“ M , iG{-M,-M + 2,-..,M-2,M} 
where L = (A/" — 1 + C/) / 2, M = (N — 1 — U)/2). The random variable U is defined as 

N—2 

U = X! a bf ' a bi+l 
1=0 

where {aij} is the signature sequence of user 1. By assuming random signature sequences, the pdf of U 
is given by 


Pu(i) 


N-l \ 

i + N-1 ) 2- 7V+1 , ie{-N + 2,-N + 3,---,N-3,N-l} 
2 / 


If 0(t) is a rectangular pulse, we have 


w k = U k S k + A k (T c - S k ) + X k T c + Y k (T c - 2 S k ) 

Let S = (Si, • • • , Sk) and $ = (</> i, • • • , 0k)- Then the conditional variance of W 1 is given by 


Var [W | S,$,M] = E 


K 


ysJP/2 ■ cos P ^ W k ■ cos (f> k ) | S, 3>, M 


fc= 2 


= y COS 2 p Y, E i W k I Sfc’ M ] ■ E [ c °s 2 fa | 0fc] 


if 


k= 2 


iC 


= ^ cos 2 P ^ [1 + cos(20fc)] • Var [W fc | S k ,M] 
k—2 


(B-6) 
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The conditional variance of W k can be computed as 


Var [W k | S k , M] = E [U 2 S 2 \ S k ] + E [A 2 {T C - S k ) 2 | S k ] + E [X 2 T 2 \ M] + E [Y 2 (T c - 2 S k ) 2 


The random variables 11 ^ and A k have variances equal to 1 . Then we have 

e [n ls 2 k | s k ] = s 2 

E [A 2 k (T c - S k ) 2 | S k \ = (T c - S k ) 2 
E [XlT 2 | M\ = T 2 (N - M - 1) 

E [Y 2 (T c - 2 S k ) 2 | S k , M] = M (T c - 2 S k ) 2 
Substituting Eqs. (B-8) through (B-ll) in Eq. (B-7) gives 

Var [W k | S k , M } = 2(2 M + l)(S 2 k - T c S k ) + NT 2 

and thus from Eq. (B-6), 


K 


Var [tT | S, M\ = | cos 2 /?E I [l + cos(20 fc )] [2(2M + l)(S 2 k - T c S k ) + NT 2 

k—2 


= -w c °s 2 p ^2 C 1 + cos ( 2 ^)] 


K 


k—2 


(2M + l)(S 2 k - T c S k ) + 


NT 2 


By averaging over (j> k , 


P 


K 


Var [W | S, M\ = — cos 2 /?E 


k—2 


(2M+l)(Sl-T c S k ) + 


NT 2 


j cos2 1 3 


K 


^(2M + l )(Sl-T c S k ) 


\_k—2 


(K - l)NPT 2 


cos 2 (3 


By averaging over S k , since E[S k — T c S k ] = —T 2 / 6 , we have 


Var [IV | M] 




PT 2 

C 

12 


K 

cos 2 (3 E( 2M+1 ) 

k—2 


I 5 fc ,M] 
(B-7) 

(B-8) 

(B-9) 

(B-10) 

(B-ll) 
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For random signature sequences, E[M ] = (N — l)/2; thus, we have 


Var [W] 




PT 2 

C 

12 


K 

COS 2 / ? 

k—2 


(K - 1)NPT 2 cos 2 /3 

~ 6 

Therefore, the variance of IT 7 is given by 


Var [W 1 ] 


(. K - 1 )NPT 2 cos 2 /? 

6 


Similarly, it can be shown that 


Var[lV Q ] 


(K - 1 )NPT 2 sin 2 (3 

6 
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